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We investigate the influence of convection on the formation of molecular spectral lines in the 
atmospheres of late-type giants. For this purpose we use the 3D hydrodynamical C0 5 B0LD and 
classical ID LHD stellar atmosphere codes and synthesize a number of fictitious lines belonging 
to a number of astrophysically relevant molecules, C2, CH, CN, NH, OH. We find that differences 
between the abundances obtained from molecular lines using the 3D and ID model atmospheres 
are generally small at [M/H]=0.0, but they quickly increase at sub-solar metallicities and may 
reach ~ — 0.9dex at [M/H]=-3.0. The 3D-1D abundance corrections show a significant depen- 
dence on the spectral line parameters, such as wavelength and excitation potential. Our compar- 
ison, therefore, points to a complex interplay between the spectral line formation and convection 
that should be properly adressed in stellar abundance analysis. 
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1. Introduction 

Late-type giant stars are important tracers of intermediate and old (>1 Gyr) stellar populations. 
They are luminous and can be detected in distant stellar systems of the Milky Way and Local Group 
galaxies and may thus provide important clues about the origins of different chemical species in 
the Universe. The latter can be achieved using accurate chemical abundances determined with 
realistic stellar model atmospheres and spectral synthesis calculations. However, suitability of the 
currently used ID model atmospheres for achieving this goal has been seriously questioned with the 
advent of the new generation 3D hydrodynamical model atmospheres. Doubts are raised by the fact 
that convection influences strongly the average temperature stratification and produces horizontal 
temperature and density fluctuations that are absent in the ID models, which may thus influence 
the physical conditions in the line forming regions (e.g., [3]). 

2. Aim and tools 

In this work we investigate the influence of convection on the formation properties of molec- 
ular lines in the atmospheres of late-type giants. We used 3D and ID stellar model atmospheres 
to perform the analysis of artificial spectral lines corresponding to a number of astrophysically im- 
portant molecules. The goal of this study was to derive quantitative constraints on the abundance 
estimates obtained from molecular lines using the 3D and ID stellar model atmospheres. We will 
further call the difference between the abundances obtained using the two kinds of atmosphere 
models a 3D-1D abundance correction (e.g., [4], [2]). 

Our analysis was based on the 3D and ID model atmospheres calculated using the C0 5 BOLD 
and LHD stellar atmosphere codes, respectively. Both stellar atmosphere codes share the same 
micro-physics (opacities, equation of state) and atmospheric parameters (Table [IJ). Using the 3D 
and ID model atmospheres we synthesized a set of fictitious spectral lines corresponding to sev- 
eral molecules: C2, CH, CN, NH, OH. Spectral line synthesis was performed with the Linfor3D 
code under the assumption of LTE (for the description of all three codes see, e.g., [1], [5]). The 
wavelengths and excitation potentials of the fictitious spectral lines were 4000, 8500, 16000 A and 
0, 2 eV, respectively. The bluest wavelength was chosen as representative of molecular lines ob- 
served in this spectral region (e.g., CH G band at ~ 4300A). The 8500 A and 16000 A wavelengths 
coincide with the maximum and minimum of the H ion bound-free opacity, respectively. 

The 3D- ID abundance corrections obtained for a given spectral line with the 3D and ID model 
atmospheres were derived using the curves of growth. Only weak spectral lines (~ <5 mA) were 
used in the analysis in order to eliminate the influence of the microturbulent velocity in the ID 
models. 

3. Results and discussion 

3.1 3D-1D abundance corrections 

The 3D-1D abundance corrections derived at four different metallicities from various molecu- 
lar lines with different wavelength, A, excitation potential of the lower level % are shown in Fig. 
The influence of convection is generally largest at lowest metallicity ([M/H] = -3.0). At a given 
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573 x 573 x 243 
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-2.0 
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Table 1: Parameters of the 3D CO BOLD models of late-type giants used in the present study. 
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Figure 1: The 3D-1D abundance corrections for different molecular species derived using fictitious lines 
with excitation potentials % = 0, 2 eV, plotted versus metallicity at three different wavelengths, 4000, 8500, 
16000 A. 



metallicity, the size of 3D-1D abundance corrections decreases with increasing excitation potential 
of a given spectral line, % (Fig. The 3D-1D abundance corrections at low metallicities differ 
significantly in their magnitude which ranges from ~ -0.1 dex in case of CH to ~ -0.9 dex in case 
of CN. 

3.2 Thermal structures in the 3D models 

The comparison of the temperature profiles of the 3D model with the corresponding temper- 
ature stratifications of the ID LHD model and the temporarily and spatially averaged 3D hydro- 
dynamical (hereafter, (3D)) model is shown in Fig. ||. At solar metallicity the ID LHD model 
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Figure 2: Temperature profiles of the 3D hydrodynamical and classical ID models of the late type giant 
plotted versus Rosseland optical depth at two different metallicities, [M/H] = 0.0 (left) and [M/H] = -3.0 
(right). The contour plot shows the density of the individual temperature profiles of the 3D model. The 
dashed line represents the temperature profile of the temporally and spatially averaged (3D) model, while 
the solid line represents the corresponding ID model. Horizontal black lines in the upper part of the plots 
show the approximate location where various molecular lines form in the 3D and ID models (solid and 
dashed lines, respectively). 



closely traces the most probable temperature profile of the 3D model. However, at [M/H] = -3.0 
and beyond log Tr oss ~ —3.5 the 3D model is cooler than the corresponding ID model (Fig. |||). 

The formation of different molecular lines occurs at different optical depths, both at solar 
and sub-solar metallicities (Fig. ^|). However, since the differences between the temperature strat- 
ifications of the averaged 3D and ID models are small at [M/H]=0.0 the corresponding 3D-1D 
abundance corrections are minor. The situation changes at lower metallicities where even small 
differences in the location of molecular line formation are associated with significant changes in 
the local temperatures, due to rapidly increasing differences between the temperature profiles of 
the averaged 3D and ID models in the outer atmosphere. This may explain why, for example, the 
3D-1D abundance corrections are smallest in the case of C2 and CH, which form deepest in the 3D 
atmosphere and are least affected by differences in the temperature profiles. 

3.3 Comparison with the work of other authors 

The only similar study where the effects of convection on spectral line formation in late-type 
giant stars were analyzed was done by [5]. The results obtained in the two studies agree in a sense 
that largest 3D-1D abundance corrections are seen in case of the most metal- poor atmosphere 
models (Fig. |3]). However, the magnitude of the abundance corrections obtained in our study is 
significantly smaller. The origin of these differences is not fully clear yet. In part, it might be 
related to the fact that the (3D)-1D temperature differences are larger in [5], particularly at lowest 
metallicites, but one should keep in mind that gravities in the two sets of models are also slightly 
different. 
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Figure 3: The 3D-1D abundance corrections for CH 4360 A, NH 3360 A, OH 3150 A lines. Dashed lines 
are results of [3], solid lines represent calculations obtained in this study. Table inside the figure gives the 
parameters of the 3D models used by [3]. 



4. Conclusions 

Convection alters significantly the thermodynamical structure of late-type giant atmospheres, 
especially in the outer layers, and its influence grows larger with decreasing metallicity. At lowest 
metallicities, the temperature in the outer layers of the averaged 3D model is significantly lower 
than in the corresponding ID model atmosphere. This leads to significant differences between the 
molecular abundances derived using the 3D and ID model atmospheres: for example, in case of 
CN they may be as large as -0.9 dex at [M/H] = -3.0. The 3D- ID corrections predicted for the 
individual molecules are also markedly different, with the largest spread at lowest metallicities. 
The strong dependence of the 3D-1D abundance correnctions on oscillator strength, wavelength, 
and excitation potential of particular spectral line clearly exposes the limitations of the ID models 
in stellar abundance work, especially when studying metal-poor stellar populations. 
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